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Abstract 

We establish empirical quantile process CLTs based on n inde- 
pendent copies of a stochastic process {Xt : t £ E} that are uni- 
form in t & E and quantile levels a G I, where / is a closed sub- 
interval of (0,1). Typically E = [0,T], or a finite product of such 
intervals. Also included are CLT's for the empirical process based on 
{Ixt<v ~ Pr(Xt < y) : t e E, y G R} that are uniform in t € E,y G R. 
The process {Xt : t e E} may be chosen from a broad collection of 
Gaussian processes, compound Poisson processes, stationary indepen- 
dent increment stable processes, and martingales. 



1 Introduction 



Let X = {X{t) : t £ E} he a stochastic process with P{X{-) G D{E)) = 
1, where S is a set and D{E) is a collection of real valued functions on 
E. Also, let C = {Cs,x- s e e M}, where Cs,x = {z & D{E): z{s) < 

^ are i.i.d. copies of the stochastic process X and 
F{t,x) := P{X{t) < x) = P{X{-) e Ct^x), then the empirical distributions 
built on C (or built on the process X) are defined by 

^ n 1 " 

Fn{t,x) = -Y,^i-oo,x](^t{t)) = -Y,I{X,eCt,.},Ct,xeC, 
1=1 i=l 

and we say X is the input process. 

The empirical processes indexed by C (or just E xM) and built from the 
process, X, are given by 

Un{t, x) := y/n{Fn{t, x) - Fit, x)) . 
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In [KKZlOj we studied the central limit theorem in this setting, that is, 
we found sufficient conditions for a pair {C,P), where P is the law of X on 
D{E), ensuring that the sequence of empirical processes {z/„(t,x) : {t,x) G 
E X M}, n > 1, converge to a centered Gaussian process, G = {Gt^x ■ (t, x) G 
E xM.} with covariance 

E{G{s,x)G{t,y)) = E{[I{Xs < x) - P(X, < x)][I{Xt < y) - P{Xt < y)]) 

This requires that the law of G on iooiE x M) (with the usual sup-norm) 
be Radon, or equivalently, (see Example 1.5.10 in jvd VW96] ) . that G has 
sample paths which are bounded and uniformly continuous on E xM. with 
respect to the psuedo-metric 

d{{s,x),{t,y)) = E{[I{X, < x) - I{Xt < y) - {P{X, < x) - P{Xt < y))f)l 

(1) 

It also requires that for every bounded, continuous F: £oo{E x M) — M, 

lim E*F(zy„) = EF(G), 

n—>-co 

where E* denotes the upper expectation (see, e.g., p. 94 in |Dud99j ). 

The quantiles and empirical quantiles are defined as the left-continuous 
inverses of F(t,x) and Fn(t,x) in the variable x, respectively: 

r^(t) = F-^{t,a) = inf{x: F{t,x) > a} (2) 

and 

Tait) = F-\t,a) = inf{x: F„(t,x) > a}. (3) 
The empirical quantile processes are defined as 

V^{F-\t,a)-F-\t,a)), 

and we also use the more compact notation 

for these processes. Since we are seeking limit theorems with non-degenerate 
Gaussian limits, it is appropriate to mention that for a e (0, 1) and t fixed, 
that is, for a one-dimensional situation, a necessary condition for the weak 
convergence of 

V^ir^it) - T^{t)) ^ (4) 

where ^ has a strictly increasing, continuous distribution, is that the dis- 
tribution function F{t,-) be differentiable at Ta{t) and F'{t,Ta{t)) > 0. 
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Hence F{t,-) is strictly increasing near Ta{t) as a function of x, and if we 
keep t fixed, but ask that holds for all a £ (0,1), then F{t,x) will 
be differentiable, with strictly positive derivative F'{t,x) on the the set 
Jt = {x :0 < F{t,x) < 1}. Moreover, by Theorem 8.21, p 168, of |Rud66j . 
if F'(t,x) is locally in Li with respect to Lebesgue measure on J^, then 
F'{t,x) is the density of F{t, •) and it is strictly positive on Jf. For many 
of the base processes we study here, = R for all t £ E, but should that 
not be the case, it can always be arranged by adding an independent ran- 
dom variable Z with strictly positive density to our base process in order to 
have a suitable input process. In p articular, the reader should consider a 
base process as one which, after possibly some modification, will be a suit- 
able input process. At first glance perhaps this may seem like a convenient 
shortcut, but we know from [KKZIO] that when E = [0, T] and the base 
process is a fractional Brownian motion starting at zero when t = 0, then 
the empirical CLT over C fails, but by adding Z as indicated above it will 
hold. In these cases adding Z is just starting the process with Z. Hence a 
typical assumption throughout sections 2, 3, and 4 will be that the distribu- 
tions F{t, •) are continuous and strictly increasing on R, but in section 5 we 
show how to remove this assumption for fractional Brownian motions and 
symmetric stable processes with stationary independent increments when 
E = [0, T], and the processes start at zero at t = 0. 

In section [2] we extend a result of Vervaat |Ver72| on the relation between 
empirical and quantile processes. Because of this relationship, when such a 
theorem is applied to empirical quantile processes built on a stochastic pro- 
cess, of necessity, the hypotheses includes the CLT for the empirical process 
as well as the existence of and conditions on the densities of F{t, x). Hence, 
in order to prove empirical quantile CLT's for stochastic processes we will 
need the associated empirical CLT's. Some of these CLT's follow from the 
results in [KKZIJD], but in section [3] we also show that several other classes of 
processes fall within the scope of those results. This includes the cases when 
the base process is a strictly stable process with stationary and independent 
increments, certain martingales, and even other independent increment pro- 
cesses. Section 4 then turns to the task of obtaining the empirical quantile 
process CLT's for t hese examples, and as mentioned above, section 5 looks 
at empirical quantile results for some important examples where one can 
get around the difficulty imposed when the input process starts zero at time 
zero. The results of section 5 were motivated by the CLT for the sample 
median of independent Brownian motions with value at 0, a result of J. 
Swanson ( [Swa07] ) , and we extend that result by proving a CLT for such em- 
pirical quantiles uniform, not only in the time parameter, but also uniform 
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in all quantiles. Moreover, we do this for symmetric independent increment 
stable processes and fractional Brownian motions. 

2 Vervaat's Approach 

Throughout we assume the notation of section one. In particular, in this 
section we are assuming that for all t £ E, F{t, x) is strictly increasing and 
continuous in x S M. Our goal in this section is to prove an analogue of 
Vervaat's Lemma 1 in |Ver72] . We follow Vervaat's idea of using an almost 
sure version of the empirical CLT. 

Notation 1. For a function / : S — > M we use the notation /* to denote a 
measurable cover function (see Lemma 1.2.1 [vdVWQG] ). 

An important result regarding weak convergence or convergence in law, 
in a general context, is that it also has a form allowing almost sure conver- 
gence. Such results have a long history, and here we use Theorem 3.5.1 in 
|Dud99] . which is slightly more general. 

Theorem 1 ( jDud99j ). Let {D, doo) be a metric space, (il. A, Q) be a proba- 
bility space and fn- — )■ -D for each n = 0, 1, • • • . Suppose /o has separable 
range, Dq, and is measurable with respect to the Borel sigma algebra on 
Do- Then {/„: n > 1} converges weakly, or in law, to /o iff there exists a 
probability space {0,,J^, P) and perfect measurable functions Qn from (fi, 
to (O, ^) for n = 0, 1, • • • , such that 

Pog-^=QonA (5) 

for each n, and 

dl^ifn o gn, fo o go) ^ 0- (6) 

a.s 

where d*^{fn ° 9n, fo° do) denotes the measurable cover function for doo(/n ° 
9n, fo ° go) and the a.s. convergence is with respect to P. 

In our setting the metric space D is loo{E x M), with distance doo the 
usual sup-norm there, and the probability space (O, A, Q) supports the i.i.d. 
sequence {Xj : j > 1} and the Gaussian process G. Then, for u; G fi, n > 1, 
the /„ of Dudley's result is our 

Ui^) = MFni; ■){oj) - F{; •)) €(oc{Ex R), (7) 

and {fn- n > 1} converges in law to 

fo{io) = G(;-){uj)(^ioo{ExR). (8) 
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That is, we are assuming the empirical CLT over C, and therefore Theorem 
[1] imphes there is a suitable probability space {^},J-',P) and a set ili C $7 
with P{Qi) = 1 such that for all a) S ili, 

||(/nO5n)H-(/0O50)Hir = ( SUp | (/„05„) (ci) - (/oO^q) M | )* ^ 0. (9) 

teE,xeR 

Hence, if 

Fn{t,x){u}) = Fn{t,x){gn{uj)) and G{t,x){uj) = G{t,x){go{u;)), 

then on r^i we have the empirical distribution functions {Fn : n > 1} satis- 
fying 

\\V^{Fn-F)-G\\* = { sup \V^{Fn{t,x){u:)-F{t,x))-G{t,x){Lj)\y ^0. 
teE,xeR 

(10) 

Remark 1. The functions Fn are still distribution functions as functions of 
X, and on Q we have ^/n{Fn — F) — G G ioo{E x E.). In addition, since 
the functions {gn : n > 0} are perfect and ([5]) holds, it follows for every 
bounded, real valued function h on ioo{E x M), and n > 1, that 

E1[MV^(F„-F))] =E^[/i(V^(F„-F))] andEp[MG)] =EQ[/i(G)]. (11) 

Since we are assuming {-y/n(F„ — F): n > 1} converges weakly to the 
Gaussian limit G, and G has separable support in ioa{E x M), then (jllh 
immediately implies {^/n{Fn — F): n > 1} also converges weakly to G. 

The generalized inverse of Fn{t, ■) in the second variable is given by 

?^{t) = F-^{t,a) = -m.i{x: F„(t, x) > a}, t G F, a G (0, 1), n > 1, (12) 

and as before for each t €z E,a £ (0, 1), the inverse function 

Ta{t) = F^^{t, a) = inf{x : F{t, x) > a}. (13) 

Of course, since we are assuming F(t, x) is strictly increasing, this is a 
classical inverse function, and to emphasize that the inverse is only on the 
second variable we also will write F~l and Ff^ for these inverses. Then, 

for each t G E we have F^l{-) : [0, 1] M and since Ft is assumed con- 
tinuous we have Ff'^i-) : (0, 1) ^ M. It is also useful to define ^^"^(0) = 
-oo,F4(-oo) = F„,j(-oo) = 0,Ft-^(l)_= oo,Ft(+oo) = F„,t(+oo) = 1, and 
Gt(-oo) = Gt{+oo) = 0. We also set 1 = M U {-oo} U {+oo}. 

To use (jlOp we will need the function Fn^t ° Ff^ and its inverse, which 
is determined in the next lemma. 



5 



Lemma 1. For each t G E 

{Fn,toF^^)-^ = FtoF-} 
where the inverses are defined as in p2|) and ([13]). 



(14) 



Proof. For each t G E,a G [0,1], we have, since we are assuming Ft{-) is 
strictly increasing and continuous, that 

{Fn,t oFi^)-\a) = inf{/3: Fn,t o F,-i(/3) > a} 
= mf{Ft{x): Fn,t{x)>a} 
= Ft{mf{x: Fn,t{x) > a}) 
= {FtoF-l){a). □ 

The next lemma is our modification of Lemma 1 in [Ver72| applicable to 
the present situation. 

Lemma 2. Let a„ — )• 0, assume that uniformly in t G E, Gt{Ff^ (a)) is a 
uniformly continuous function of a G (0, 1), and 



sup 



Fn,t{x) - Ft{x) 



Gt{x) 



(15) 



as n tends to infinity. Then, setting /t(a) = a for t £ E,a £ [0, 1], we have 



sup 

, ieB,Qe[o,i] 



{Fn,toFf'){a)-It{a) 



iGtoFi^){a 



and 



sup 

[0,1] 



{FtoF-}){u)-It{u) 



{GtoF-\u) 



0. 



(16) 



(17) 



Proof. Since we are assuming for each t £ E that F{t, •) is strictly increasing 
and continuous on M, it follows that {Ff (a) : a G (0, 1)} = M. Therefore, 
if one restricts a in (|16p to be in (0, 1), then (jlOp follows immediately from 
(jl5p . To obtain (|16p for a = and a = 1, then follows from the conventions 
we made prior to the statement of the lemma involving ±00. 

To show ()16p implies (|17p we define for each t £ E the completed graph 
of {Fn,t o I^t''^){-) on [0, 1] to be given by 



Tn,t = {{a,u): a £ [0, 1], (F„,t o Ff^){a - 0) < u < {Fn,t o Ff^){a + 0)}. 
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Here {Fn,t°Ff^){a±0) denotes the left and right hand limits of {Fn,t°Ff^){-) 
when Q G (0, 1), and are given through the conventions above when a = 
or 1, i.e. we understand the left and right hand limits at zero to both be 
zero, and the left and right hand limits at one to both be one. Now (jl6p 
implies that 

Jim^(sup{|^^-(GtoF-i)(a)| : {t, a) G E x [0,1], {a,u) gT^^^})* =0, 
which also implies that 

^lini^(sup{|^ + (GtoFr')(a)| : {t,a) £ E x [0,l],(a,n) G r„,t})* = 0. 

(18) 

For each t G E we set 

r-] = {{u, a): ue [0, 1], {Ft o F~J){u - 0) < a < {Ft o F-}){u + 0)}, 

where one can check that the left hand limit of Ft o F^l{-) at zero is zero 
and we take the right hand limit at one to be one. 
Then one can check that 

{a,u) G r^^t if and only if {u,a) G r~j. (19) 

Moreover, ()17p is implied by 



Jim (sup{|^ + {Gt o F^^){u)\ : {t,u) e E x [0, 1], (n,a) G r;]})* = 0, 

(20) 

and (jlSp and (jl9p implies (j20p provided we show 

lmi{sup{[{GtoFi'){u)-{GtoFf'){a)[: t € E,{u,a) G T-},})* = 0. (21) 

Since we are assuming the empirical CLT holds over C with Gaussian 
limit process {G{t, x) : {t,x) £ E x R}, it follows that G has a version which 
is sample uniformly continuous on E x M with respect to its L2-distance, 
d{-,-), given in ([1]). This is a consequence of the addendum to Theorem 
1.5.7 of jvdVW96] . p. 37. When referring to G we will mean this version. 
By the total boundedness of the distance, the associated space of uniformly 
continuous functions is separable in the uniform topology. This space is a 
closed subspace of i^o {E xM), which then implies that this version of G is 
measurable with respect to the Borel sets of ioo{E x R). Using the definition 
of G following Q, and ([5]) with n = 0, we have the laws of G and G are 
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equal on ioo{E x M). In particular, G is also measurable to the Borel sets 
of ioo{E X M) with separable support there, it has the same covariance G 
and its L2-distance is d, and it is sample continuous on (E x M, d) with 
P-probability one. 

Now for each t £ E and a, /? G (0, 1) we have 

d{{t,Ff\a), it,Ff\p)) = \a-(3\-\a- . (22) 

Thus for each t £ E we have d{{t,Ff^{a), {t,Ff^{f3)) as q, /3 ^ or 
a, /3 — )■ 1. We also have 

d{{t,Ff^{a),{t,Ff^{0)) = 0-0^ < |a-0|, 
d{{t,Ff^{a),{t,Ff^{l)) = 0-0^ < |a-l|, 

and 

d{{t,Fi\0),{t,Fi\l)) = 0, 

and hence the uniform continuity of G along with G{t, Ff^{0)) = G{t, Ff^{l)) 
= implies that uniformly in t G we have Gt{F^^ (a)) is uniformly contin- 
uous in a € [0, 1] with probability one. Moreover, the process {G{t, x) : (t, x) 
€ E X M} has separable support in £^(^E X M), and hence the upper cover 
used in (f2T]l is unnecessary as the function there is measurable. We also 
have with P-probability one that 

sup \Gt{Fi\a))\ < oo, (23) 
te-B,oe[o,i] 

and hence converging to zero, and (fTSj) implies 

lim sup{|ii — a\: t £ E, (a, u) £ Tnt} = 0. (24) 

n— >oo ' 

Therefore, (|2ip follows from (|19p and that uniformly in t G ii^ we have 
Gt o Ff (a) uniformly continuous in q G -B x [0, 1]. Therefore, (I2ip holds, 
and this implies ()17p . so the lemma is proven. □ 

2.1 Applying Lemma [2] to an Empirical CLT 

Assuming the empirical CLT holds over C, the conclusions of Lemma [2] 
hold with a„ = -7=, and we have proved the following lemma. 
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Lemma 3. For sdl t G E assume the distribution function F{t, x) is strictly 
increasing and continuous in x G M, and that the CLT holds on C with limit 
{G{t,x): {t,x) £ E X M}. Then, with P-probability one, we have 

sup \V^[{Ft o F-}){a) - Itia)] + (Gt o ^ 0. (25) 

ie£,ae[o,i] ' 

Up to this point we have only assumed that the distribution functions 
{Ft{-): t G E} are continuous and strictly increasing on M, and that the 
empirical processes satisfy the CLT over C. Now we add the assumptions 
that these distribution functions have densities {f{t, ■): t G E} such that 

limsup sup \f{t,u)- f{t,v)\=0, (26) 

o-^O teE \u~v\<5 

and for every closed interval I in (0, 1) there is an 9{I) > such that 

inf f{t,x) = cj0(n>O. (27) 

Lemma 4. Assume for all t G that the distribution functions F[t, x) are 
strictly increasing and continuous, and that their densities f{t,-) satisfy (I26p 
and (|27p . If the CLT holds on C, then for every closed subinterval / of (0, 1) 

lim[ sup |f^(t) -T„(t)|]* =0. (28) 

in P probability. 

Proof. Since we are assuming (j27p . fix / a closed subinterval of (0, 1), and 
take < e < 6'(/). Let 

An = {[ sup \?^{t)-T^{t)\r >e}, (29) 
teE,aei 

and 

5„ = {[ sup \Fr,{t,x)-F{t,x)\]* >6}, (30) 
teE,xeM. 

where < 6 < 5{e) < Then, since we have the CLT over C with 

respect to P, Lemma 2.10.14 on page 194 of jvdVW96] implies there exists 
ns < oo such that n > ns implies 

PiBn) < e. (31) 
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In addition, (|26p and (|27p imply we also have 



and 



sup F{t, Ta{t) -e) <a-6, (32) 



inf rc,(t) + e) > a + 5. (33) 

a£l,t£E 



That is, (132]) holds by (I2ZD since 

r-TQ (t) 



sup F(t,Ta{t) — e) < a — inf / f{t,x)dx, 



and if (5 = 5(e) < ^^^,0 < e < 9{I), we then have 



/■Ta (t) 

inf / f{t,x)dx > ecj QU) > 6. 



r-Ta (t) 
'Tc{t)-e 



Thus on B^, for all t S a G /, 

F(t,T^(t)) >Fn(t,T^(t))-5>a-<5, 

where the second inequality follows by definition of T^{t). Combined with 
P2p . on this implies that for all t G -E, all a G / 

T^{t)>T^{t)-e. (34) 

Similarly, (|33l) holds by (|27l) since 

fT-a(t)+e 

f Fft.T„m + e) =Q+ inf / 
and if (5 = 5(e) < ^^^,0 < e < 0(1), we then have 

(■Ta{t) + <i 

Thus for X <T2{t) and alH G i?, a G /, on we have 

F{t,x) < Fn{t,x) + 6< a + 6 < F{t,Ta{t) + e), 

where the first inequality follows from the definition of B!j^, the second by 
definition of r^(t) and that x <T^{t), and the third by ([33|) . Thus Ta{t)+e > 
X for all X <T2{t), and on B^ we have 

T„(t) + e>T2(t) 



inf F(t,Ta(t) + e) = a + inf / f (t, x)(ix > a + 5 

f^,0 < e < 0{I), we then have 
inf / f(t,x)dx>ecTg(T\>5. 
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for all t £ E,a G I. Combining this with (j34p . on we have for all t G E, 
all Q G /, that 

Ta{t)-e<?^{t)<T^{t)+e. (35) 

Hence on 

[ sup K{t)-T^{t)\r <e, 

and on Bn it is certainly bounded by oo. Since B^ is measurable, we thus 
have for n > ns that 

P([ sup |f^(t)-T„(Or >e) <P(5n) <e. 

teE,aei 

Since e > can be taken arbitrarily small, letting n — t- oo implies (1280 . Thus 
the lemma is proven. □ 

Theorem 2. Assume for all t G E that the distribution functions F{t,x) 
are strictly increasing, their densities f{t,-) satisfy (j26p and (j27p . and the 
CLT holds on C with limit {G{t,x): {t,x) £ E xR}. Then, for / a closed 
subinterval of (0, 1) we have 

( sup \Mr:{t)-r^{t))f{t,T^{t)) + G{t,Ta{t)\)* ^0 (36) 

in P-probability, and therefore the quantile processes {\/n(^(0 ~ '^aifS) 
f{t, Ta{t)) : n > 1} satisfy the CLT in too{E x /) with Gaussian limit process 
{G{t, Ta{t)) : (t, a) S X /}. Moreover, the quantile processes {■\/niT2{t) — 
Ta{t)) : n > 1} also satisfy the CLT in iao (E x /) with Gaussian limit process 

{ffei^ (t,a)Gi?x/}. 

Proof. Applying Theorem 3.6.1 of |Dud99j . the first CLT asserted follows 
immediately from (j36p . Hence we next turn to the proof of (j36p . 

First we observe that under the given assumptions, we have ()25p holding. 
Furthermore, since the densities are assumed continuous, 

F{t, y) - F{t, x) = fit, x){y - x) + R{t, x, y){y - x), (37) 

where R{t, x, y) = f{t,(^(t)) — f(t, x), and ^(t) between x and y is determined 
by the mean value theorem applied to F{t, •). Of course, R{t, •, •) depends 
on F(t, •), but we suppress that, and simply note that since ^(t) is between 
X and y, 

\R{t,x,y)\< sup \f{t,u)-f{t,x)\. (38) 

uG[x,y]U[y,x] 
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Therefore, for M > we have 

P{[ sup \^{?^(t)-T^ml^^^^^]*>M)<an{M) + bn, 

where 

a„(M) = P(^„,,in^„,,2), 

A,,i = {[ sup |V^(?^(t)-r,(t))|(/(t,r,(t)) + i?(t,f^(t),r„(t)))]* >M}, 

A„,2 = {[ sup |i?(t,f^(t),r„(t))r<^}, 

bn = P{[ sup |i?(t,f^(t),T„(t))|]* > 

and Cj 0(^j-^ > is given as in ([27|) . Thus 

P([ sup |V^(?^(t)-r„(t))|:^^^i^^]*>M)<P(^„,i)+6„, 

and by (|37|) we also have 

A„,i = {[ sup \V^{Ft{?^{t)) - Ftir^mir > M}, 

t£E,a£l 



which imphes 

P([ sup \V^{9:{t) - r„(t))|:^^^l^]* > M) (39) 

t&E,a&I ^ 

< P{[ sup |^/^(Fi(f^(t))-Fi(r„(t)))|r >M) + 6„. 

Since It{a) = Ft{Ff^{a)),a G (0, l),t E ^, and / C (0, 1) 

[ sup \V^{F,{?^{t)) - Ftir^mir (40) 

t&E,a&I 

< [ sup \MiFtoF;;}){a)-It{a)] + {GtoFr\a))\r 

teE,ae[o,i] 

+ [ sup \GtoF,-\a))\r, 

t&E,a&I 

and since the process {G{t, x) : t £ E,x M} is sample continuous on x M 
in the semi- metric d given in ([TJ with Radon support in ioo{E x M) we also 
have 

[ sup \GtoFf\a)\]* = sup \GtoFi\a)\. (41) 
teE.aei teE,aei 
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Therefore, for every e > and all n > 1, by combining (j25p . ()40p and (j4ip 
we have an M = M(e) sufficiently large that 

P{[ sup \V^{F,{T^{t)) - F,{r^m\r >M)<e. (42) 

t&E,a&I 

We now turn to showing that 

[ sup |V^(f^(i)-r„(t))|]* (43) 

is bounded in P-probability. That is, let 

\{t,5)= sup \f{t,u)- f{t,v)\. 

\u—v\<5 

Then 

\R{t,x,y)\ < X{t,\x-y\), 

and hence by (j26p for every e > there exists 6 > such that |x — y| < 6 
implies 

snp \R{t,x,y)\ < e. 

Therefore, for every e G (0, -^^|^) there exists 6 = 6{e) > suitably chosen 
such that 

6„<P([ sup \R{t,?:{t),T^mr>e) (44) 

t£E,a£l 

= P*i sup |i2(t,f^(t),r„(t))| >e), 

and since 

P*{ sup |i2(t,f^(t),r„(t))|>e)<P*( sup \?-{t) - r^t))] > 6) (45) 

t£E,a£l teE,aeI 

= P{[ sup |f:(t)-r„(t))|r ><5), 

Lemma m implies for every e E (0, -^^|^) that 

hm bn = 0. (46) 

Combining ([MI), (112]), and (gSD, we have dS]), i.e. [sup^g^^^g^ \Vn{^it) - 
Ta{t))\]* is bounded in P-probability. Furthermore, we then also have that 

[ sup \v^i?:it)-T^m\Rit,?:{t),T^m]* m 

t&E,a&I 
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converges in P probability to zero. 

Now, by (j25p and ()37p we have with P-probabihty one that 



Hm sup |V^(?^(t)-T«(t))[/(i,ra(t))+i2(t,r„(f),f^(t))]+G(t,ra(t)| 

n^'^\t(iE,a&I 

(48) 

and since 



[ sup |V^(f^(t) -r„(t))/(t,r„(t)) + G(t,T„(t)|]* <u„ + ^;„, 



where 



un < [ sup \v^{?:^{t)-Tamf{t,T^{t))+R{t,TUt),9:m+G{t,TUt)\y 

t&E,aeI 

and 

vn<[ sup |^/^(f:(^)-r„(^))i^(^,r„(^),f:(^))|]^ 

we have by combining (j47p and (j48p that 

[ sup \^(9^{t)-r^it))fit,T^it)) + Git,T^{t)\]* ^0 



in P probabihty. Hence ()36p is proven. 

To finish the proof it remains to check that the quantile processes 

{V^{T2{t)-T^{t)): n>l}. 

also satisfy the CLT in £00 {E x /) with Gaussian limit process { : {t,a 

E X I}. Since (136p holds, and by (1271) we have the non-random quantity 

1 

sup -F7T TTTV < 

t&E,aeI j[t,Ta[t)) 

we thus have 

( sup |v^(r^(t)-r.(t))+ ^//'^f/^ )%0 (49) 

^teE,aeI flt.Tr.it]) / 



in P-probability. The CLT then follows from Theorem 3.6.1 of |Dud99j . and 
that the Gaussian process G is symmetric. Hence the theorem is proven. □ 

The next result shows that the conclusions of Theorem [2] also hold for 
the relevant processes as defined on the original probability space {i},A, Q). 
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Corollary 1. Assume for allt & E that the distribution functions F{t, x) are 
strictly increasing, their densities /(t, •) satisfy ([26l) and (f27|) . and the CLT 
holds on C with limit {G{t, x) : {t, x) G £^xM}. Then, for / a closed subinter- 
val of (0, 1), the quantile processes {■\/n{T^{t)—Ta{t))f{t, Ta{t)) : n > 1} sat- 
isfy the CLT in £oo{E x I) with Gaussian limit process {G{t, Ta{t)) : {t, a) G 
E X /}. Moreover, the quantile processes {\^{T^{t) — Ta{t)): n > 1} also 
satisfy the CLT in ^oo {E x /) with Gaussian limit process { : (t, a) G 



Proof. Recall the notation established at the start of this section in connec- 
tion with the statement of Theorem [H and the perfect mappings g„ : Q — )• 
such that Q = P o g~^. In particular, equations ^ to ([13|) are relevant. 
For ui, • • • , Un G D{E) and n > 1, t G £', a G (0, 1) define 



where Ta{t) = {Ft) ^(a). Hence setting 

rn{t,a,uj) = kn{Xi,--- ,Xn,t,a){uj) = k„{Xi{-,uj), ■ ■ ■ , X„(-, a;), t, a), 
we then have 



E X /}. 



n 




-{Ft)-\a)]f{t,T^{t)). 



V^[iFn,t)-Ha){u;) 



Ta{t)]f{t,Ta{t)) 



rn{t,a,uj) 



(51) 



and 



Ta{t)]f{t,Ta{t)) 



(52) 



Therefore, for w G 




h{,/^[{Fn,t) \a){u)-Ta{t)]f{t,Ta{t))) = (/i o r„ (t, a, •) o 5„) (w) , 



and hence the upper integrals 
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[ihornit,a,-)ogn)noj)dPiu;) (53) 
{[horn{t,a,-)]* o gn){io)dP{Lo), 
where the last equahty holds since gn is perfect. Now 



{[horn{t,a,-)]* o gn){u})dP{Lv) = / [ho rnit,a,Lv)]*dQ{uj) 



n 



n 



{h o rn){t, a,u>)dQ{uj), 

and therefore by ()5ip and (j53p . for all h bounded on loo{E x /), 

/i(V^[(i^n,t)"'(«)((i) - r„(t)]/(i,r,(t)))dP(cD) 

h{^[{F^^)-\a){u) - T^{t)]f{t,T^{t)))dQ{u). (54) 

Now the equality in ()54p implies that the quantile processes 

{V^[{Fn,t)'Ha){Co) - T^{t)]f{t,Ta{t)) : n > l,t ^ E , a ^ 1} 

satisfy the CLT in ioo{E x /) if and only if 

{V^[{Fn,t)-\a){u;) - Ta{t)]f{t,Ta{t)) : n > l,t £ E , a ^ 1} 

satisfy the CLT there, and they have the same Gaussian limit, namely 

{G{t,Ta{t)) :teE,aeI}. 

A similar argument implies the quantile processes 

{V^[{Fn,ty\a){u;) - Ta{t)] ■.n>l,teE,a£l} 

satisfy the CLT in iooiE x /) if and only if 

{V^[{Fn,tyHa){uj) - Ta{t)] ■.n>l,t€E,a€l} 

satisfy the CLT there, and they have the same Gaussian limit. Since The- 
orem [2] implies the Gaussian limit of 

{V^[(F„,t)-i(a)(w) -r„(t)] ■.n>l,teE,aeI} 
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is given by 



f{t,ra{t)) 



: t e E,a e 1} 



which has the same Radon law on £, 



■oo 



{E X /) as 



{ 



Git,Ta{t)) 

f{t,ra{t)) 



te E,ae 1} 



the corollary is proven. 



□ 



3 The Empirical CLT over C 



In order to prove the empirical quantile CLT of Theorem [2l and its corol- 
lary, we assumed the empirical CLT over C holds. Empirical results of this 
type were established in [KKZlOj for fractional Brownian motions and the 
Brownian sheet as long as these processes were not fixed to be zero at some 
point, and later in the paper we will use these facts to establish the empirical 
quantile CLT with those processes as the base process. The purpose of this 
section is to broaden the class of base processes to which [KKZIO] applies, 
and that then will also be potential applications for our quantile process 
results. In particular, in this section we show how [KKZlOj yields the em- 
pirical CLT for a broader class of Gaussian processes, all compound Poisson 
processes, and also many stationary independent increment processes and 
martingales. In particular, these results apply to all symmetric stable pro- 
cesses, and below we will show that under certain ci rcumstances empirical 
quantile CLT's also hold for such processes. 

In the typical empirical process result over C that we establish, the input 
process {X(t) : t £ E} is given in terms of a base process {Y{t) : t G E} 
where P{Y{0) = 0) = 1 when E = [0,T], and X{t) = Y{t) + Z,t e E. The 
random variable Z is assumed independent of the base process, and has a 
density which is uniformly bounded on M, or in La(M) for some a G (1, oo). 
The use of Z allows to say the densities of each X{t) have a uniform property 
provided the density of Z has that property, and is an efficient way to 
do this. More important, however, is that in many classical examples the 
base process {Y{t) : t G E} with P(Y(0) = 0) = 1 fails the empirical 
CLT over C, yet the input process X{t) = Y{t) + Z,t £ E, satisfies it. 
Examples of this type include fractional Brownian motions on E = [0,T], 
the d-parameter Brownian sheet on E = [0,T]'^, and also strictly stable 
processes with stationary independent increments. This was pointed out for 
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fractional Brownian motions in [KKZlOj . and we will say more about the 
other examples at appropriate points of this section. 

3.1 Additional Gaussian process empirical CLT's over C 

Throughout this subsection we assume E is a compact subset of the d-fold 
product of [0,r], which we denote by [0,T]'^, and that {Xt : t G E} is 
a centered Gaussian process whose L2-distance dx is such that for some 

ki < oo, s,t £ E, 

dx{s,t) = [E{{Xt-Xsf)]"^ <he^s,t), (55) 

where e(s,t) is the usual L2-distance on M"^ and < 7 < 1 . Furthermore, 
applying Theorem 6.11 on page 144 of [AGSOj we have from (j55p that {Xt : 
t G E} has a sample continuous version {Xt : t G E} such that for s,t G E 

\Xt-X,\<Te'{s,t), (56) 

where F < 00 with probability one, and < r < 7. Hence, without loss of 
generality, we may also assume throughout the sub-section that {Xt : t G E} 
is sample continuous with (|56|) holding. 

Proposition 1. Let ii^ be a compact subset of [0,r]'^, and assume {Xt : t G 
E} is a sample continuous centered Gaussian process such that ([55]) holds 

and for all x, y G M, A;2 < 00, and some /3 G (0, 1] 

sup\Ft{x)-Ft{y)\<k2\x-y\^. (57) 

t&E 

Then, the empirical CLT built from the process {Xt : t G E} holds over C. 
Moreover, if {Yt : t £ E} is a sample continuous centered Gaussian process 
such that ()55p holds, and Z is a random variable independent of {Yt : t G E} 
whose density is uniformly bounded on M, or in Lp{M) for some p G (l,oo), 
then the empirical CLT based on the process {Xt : t G E} holds over C, 
where Xt = Yt + Z,t £ E. 

Proof. First we assume {Xt : t G E} is a sample continuous centered 
Gaussian process such that (i55]l . (l56]l . and (|57l) hold. Then, applying the 
Fernique-Landau-Shepp result we have exponential decay of the tail prob- 
ability of F in (I56p . and hence assumptions (I ) and (II) of Theorem 5 in 
[KKZIO] hold. If dSSl) and §6^ hold for {Yt : t £ E} and Xt = Yt + Z, 
where the density of Z is uniformly bounded density or in Lp as indicated, 
then standard convolution formulas imply (I57p holds for {Xt : t G E}. In 
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particular, if the density is Z is assumed to be uniformly bounded, then (j57p 
holds with /? = 1, and if it is in Lp(M), then /3 = 1 — 1/p suffices. Therefore, 
under either assumption on the density of Z, we have assumptions (I ) and 
(II) of Theorem 5 in [KKZlOj holding for {Xt'.t^E}. 

Therefore, the conclusions of the proposition hold in either situation 
provided we verify condition (III) of that theorem. That is, from Remark 
8 in [ KKZIO] we need to verify there exists a centered Gaussian process 
{H°^{t) : t £ E} with L2-distance pa{s,t), which is sample bounded and 
uniformly continuous on {E,pa), and for some a G (0, ^) we have 

ieis,t)y'' <Pais,t),s,tGE. (58) 

To verify (j58p we first assume d = 1, and hence that S is a compact 
subset of [0,T]. The Gaussian process {H°'{t) : t G E} is then defined to 
be a centered, sample continuous, a6 fractional Brownian motion on [0, T] 
with L2-distance Pa{s,t) = \s — G [0,7"] and 9 G (0, r) sufficiently 

small that ([581) holds. 

If d > 1 and a G (0, ^) is fixed, then £' is a compact subset of [OjT]"^, 
and the Gaussian process {H"{t) : t G E} is defined to be the sum 

d 

H^it) = Hf{tj),t = ih, ■■■ [0, r]^ (59) 

where the are centered, independent sample continuous, a9 fractional 
Brownian motions on [0,T] such that for tj > and 6 G (0,r) 

E((/7°(tj))2) = (T V l)2™tf °. (60) 
Hence the L2-distance for {-/?"(*) : t G [0, T]'^} is 

d 

pUs,t) = (Y^iTVlfnt, - sjf "')^s,t G [0,r]^ (61) 
i=i 

and with e(s,t) the Euclidean distance on R'^, we have 

e{s,tY- = (x: \t, - s,\r-^ < iizi^^r^iT V imK m 

where the inequality holds since < ro < 1. Since 6 G (0,r) we therefore 
have 

e(s,t)™ = (X: It, - s,\r-^ < (X:(%^)''^^7^ V l)2^-)i (63) 
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Therefore, by combining ([6T]) . and ([63|) we have ([58]). and the proof is com- 
plete. □ 



Our first apphcation of Proposition [T] is to fractional Brownian motions. 
This result was obtained in |KKZ10j , but we include it here as its proof is an 
immediate application of this proposition, and an empirical quantile CLT 
result will also be obtained for these processes later in the paper. 

Corollary 2. Let E = [0, T], and assume {Yt: t £ E]} is a centered sample 
continuous 7-fractional Brownian motion for < 7 < 1 such that 10 = 
with probability one and E(y/) = t^T for t e E. Set Xt = Yt + Z, where Z is 
independent of {Yt: t £ E}, and assume Z has a density that is uniformly 
bounded on M or is in Lp{M) for some p S (1, 00). Then, the empirical CLT 
holds over C. 

Proof. The L2-distance for {Xt : t G E} is dx{s, t) = \s — tp, and hence ([551) 
holds with k — 1 = 1. Also, ()56p holds with < r < 7, and the assumptions 
on the density of Z then imply (i57|l . Therefore, Proposition [T] applies to 
complete the proof. □ 

Our next application of Proposition [T] is to the d-dimensional Brownian 
sheet. A result for d = 2 appeared in [KKZIO] . but once we have Proposi- 
tion [T] in hand, the d-dimensional case follows easily. 

Corollary 3. Let E = [0,T]'^ for d > 2, and assume {Yt: t £ E} is a 
centered sample continuous Brownian sheet with covariance function 

d 

E{YsYt) = '[l{sjAtj), s = {sir-- , Sd),t = {h, ■ ■ ■ ,td) e E. (64) 
i=i 

For t £ E, let Xt = Yt + Z, where Z is independent of {Yt: t £ E}, and 
assume Z has a density that is uniformly bounded on R or is in Lp(M) for 
some p £ (l,oo). Then, the empirical CLT based on the process {Xt : t £ 
-E} holds over C = {Ct^x'- {t^x) £ E x M}, where in this setting Ct^x = 
{z £ D{E): z{t) < x}, and D[E) denotes the continuous functions on E. 
Moreover, the empirical CLT over C fails for the base process {Y{t) : t £ E}. 

Proof. First we observe that if < Sj < tj < T for j = 1 , • • • ,d, then for 
d > 1 

d d d 

Yl*j-ll'j^^'^'J2\*^-'j\- (65) 
j=i j=i j=i 
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This elementary fact is obvious for d = 1 with = 1, and for d > 2 
it follows by an easy induction argument. Moreover, the L2-distance for 
{Xt : t G E} satisfies 

d d d 

d\{s, t) = \]\tj + J[sj-2 \[{sj A t,)l 
j=i j=i j=i 

and hence 

d d d 

d\{s, t) < T^-HY.(^3 - tj A sj) + ^{sj - tj A s,)] = T^-i \t, - s,\. 
j=i j=i j=i 

Therefore, 

which implies dx{s,t) < d^e^i [s,t), and hence ([55]) holds. Either as- 



sumption for the density of Z implies (I57p for a suitable /?, and thus Propo- 
sition [1] applies to show the CLT over C holds for {X{t) : t £ [0,T]'^} holds. 

To see why this CLT fails for the base process {Y{t) : t E [0,T]'^}, 
observe that the process W{r) = y(rd (1, • • • , 1)), < r < T'^, is a Brownian 
motion with P{W{0) = 0) = 1. Thus by Lemma 5 of [KKZlOj we have 
{Y{t) : t G [0,r]'^} fails the CLT over the the class of sets Ci = {Cr,x ■ < 
r < T'^,x e M}, where Cr,x = {z e D{E) : z{ri{l,--- ,1)) < x}. Since 
Ci ^ C, it follows from that the CLT for Y over C must also fail. □ 



3.2 Compound Poisson process empirical CLT's over C 

Here we examine the empirical CLT over C when our base process is an 
arbitrary compound Poisson process. This will be done in the next propo- 
sition by applying Theorem 3 of |KKZ10| . We will see from its proof that 
the Gaussian process needed for this application can be taken to be a sam- 
ple continuous Brownian motion, and the space of functions D{E), when 
E = [0,T], is the standard -D-space of functions on [0,T] which are right 
continuous on[0, T) and have left limits on (0, T]. These examples are some- 
what surprising since the sample paths of the base process {Y(t): t G [0,T]} 
have jumps, while those of significance in [KKZIO] and the previous subsec- 
tion were all sample path continuous. 
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To define the base process in these examples we let {A^(t) : < t < 
00} be a Poisson process with parameter A G (0,oo), and jump times 
Ti,T2,---. As usual we assume P{N{0) =0) = 1, and that the sample 
paths {N{t) : < t < 00} are right continuous and nondecreasing. Also, 
let {y^ : A; > 1} be i.i.d. real-valued random variables, independent of 
{N(t) : < t < 00}, and without mass at zero. Then, Y{t) is defined to be 
zero on [0, n), Yi on [ti,T2), and Yi -\ \-Yk on [Tk,Tk+i) for A; > 1. 

Proposition 2. The empirical process built from i.i.d. copies of the com- 
pound Poisson process {Y{t): t £ E} with parameter A € (0,oo) and 
E = [0,T] satisfies the CLT over C. 

Proof. The proof follows by applying Theorem 3 of jKKZlOj . This is ac- 
complished by showing {Y{t) : t £ E} satisfies the L condition of |KKZ10| 
when the Gaussian process involved is Brownian motion and the p distance 
is a multiple of standard Euclidean distance on [0, T]. Since the distribution 
function ofY{t) is not necessarily continuous, the L-condition involves distri- 
butional transforms of the the distribution functions F{t,x) = P{Y{t) < x) 
denoted by Ft{x). They are defined for t £ E,x £M as 

Ft{x) = F{t,x-) + V{F{t,x) - F{t,x-)), 

where ^ is a uniform random variable on [0, 1] independent of the process 
{Y{t): t£E}. 

To verify the L-condition for the Y process, let {H(t) : < t < 00} be 
a sample continuous Brownian motion with P{H{0) = 0) = 1 satisfying 

p^{s,t) =E{{H{s) - H{t) f) = 4(AVl)|t-s|. 

Then, for e > 

A = sup P( sup \Ft(Y{s)) - Ft(Ym > 

te[0,T] ^{s: p{s,t)<e} 



< sup P sup \Y{s) -Y{t)\ > 

te[0,T] ^{s: p{s,t)<e} 



sup 

te[o,T] 



1-P( sup \Y{s) -Y{t)\ = 

{s: p(s,t)<e} 



Since Y{s) - Y{t) = whenever N{s) - N{t) = 0, and for t £ [0,T] fixed 
J sup |7V(s) -7V(t)| =ol 

[{s: p{s,t)<e} J 

.2 ^ ^ // ,2 



22 



it follows that 



^l'"ll'+4(AVl)j'^^J-"U'-4(AVl)J^»J='' 

= P( sup \N{s)-N{t)\=0) 

{s: pis,t)<e} 

<P{ sup \Y{s) -Y{t)\ = 0). 

{s: p{s,t)<e} 



Now 



> exp 



Ae 



2 



2(AV1) 

and hence for < e < eo we have 



e2 



A< l-exp|--| <e'. 

Taking L suitably large we have for all e > that A < Le^, and hence 
the L-condition holds for the compound Poisson process Y, completing the 
proof of the proposition. □ 

Remark 2. Let Z he a random variable independent of both {N{t): < 
t < oo} and {1^ : A; > 1}, and define 

X{t) = Y{t) + Z,te[0,T]. 

Since the L-condition for the input process {X{t) : t G E} involves only its 
increments, and those are identical to those of the base process {Y{t) : t G 
E}, the argument above implies the L-condition also holds for {X(t): t G 
E} Therefore, the empirical process built from i.i.d. copies of X satisfies the 
CLT on C. 



3.3 Empirical process CLT's over C for other independent 
increment processes and martingales 

The processes we study here are either martingales, or stationary indepen- 
dent increment processes. There is some overlap between these examples 
and the compound Poisson processes of the previous sub-section, as such 
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processes have stationary independent increments, and could also be mar- 
tingales. However, it is easy to check that there are examples which fit into 
one and only one of the classes we study. 

Let E = [0,T], and assume {Y{t) : t > 0} is a stochastic process whose 
sample paths are right continuous , with left hand limits on [0, oo), and 
satisfying P{Y{0) = 0) = l.We say that {Y{t) : t > 0} has Lp-increments 
that are Lip-/3 on E if for some p G (0,1] and for all s,t G E there is a 
(3 £ (0, 1] and C < cxd such that 

K{\Y{t)-Y{s)\P) <C\t- s\f^. (66) 

For example, if {Y{t) : t > 0} is a strictly stable process with stationary 
independent increments and index r G (0,2] , then for r G (1,2] we have 
E(|y(t)|) = t?E(|y(l)|) and hence 

E{\Yit)-Y{s)\) = \t-s\-rE{\Y{l)\), 

which implies it has Li-increments that are Lip-^. Of course, it is also a 
martingale when r G (1, 2]. If < r < 1 , then for < p < r we have 

E(|y(t)-y(.)n = |t-.|?E(|y(i)n, 

which implies it has Lp-increments that are Lip- 2. 

li{Y{t) : t > 0} is a square integrable martingale with A(t) = E(y^(t)), t > 
0, then for < s < t the orthogonality of the increments of {Y{t) : t > 0} 
implies 

E{{Y{t)-Y{s)f) = X{t)-X{s). (67) 

Hence, if A(-) is Lip-7 on E, then (i67|) implies (f66j) with p = l,/3 = 7/2. In 
addition, if {Y{t) : t > 0} also has stationary, independent increments with 
P{Y{0) = 0) = 1 and X{t) = E(|y(t)|) < cx),t > 0, then for s,t £ E we have 

E{\Y{t) - Y{s)\) = E{\Y{\t - .1)1) < A(|t - .1). (68) 

Therefore, if A(t) < Ct^ for t G [0,<5] and some > 0, /3 G (0, 1], then it 
is easy to check that (j68p implies (j66p with p = 1 and the given (3 for all 
s,t £ E, and a possibly larger constant C. 

We also assume Z is a random variable independent of {Y{t) : t > 0} 
with density g{-) on M such that 

k = sup \g{x)\ < 00 or g G La(M) (69) 

x&M. 
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for some a £ (1, oo). Let X{t) = Z + Y{t),t > 0, and denote the distribution 
function of X{t),t > by Ft{x). Then, if g is uniformly bounded 

snv\Ft{x) - Ft{y)\ < - y|, x, y G M, (70) 

and if 5 € -La(M) we have a, k < oo 

snv\Ft{x) - Ft{y)\ < fe|x - yl^"^, x, y G M. (71) 

Proposition 3. Let E = [0, T], and assume {Y{t) : f > 0} is a stochastic 
process whose sample paths are right continuous, with left hand limits on 
[0,oo), and satisfying P(y(0) = 0) = 1. Furthermore, assume {Y{t): t G 
E} is a martingale whose Li-increments are Lip-/3 for some /3 G (0, 1] , or a 
stationary independent increments process satisfying (I66p for some p £ (0, 1) 
and /3 G (0, 1]. Let X{t) = Z + Y(t),t > 0, where Z is a random variable 
independent of {Y{t) : t > 0} and having density g{-) on R satisfying (f69]l . 
Then, the empirical process built from i.i.d. copies of {X{t) : t G E} satisfies 
the CLT over C. 

Proof. Let p{s,t) = \s — t\^ . Then, p is the L2-distance of a ^-fractional 
Brownian motion on E, and the proposition follows from Theorem 3 of 
[KKZIO] provided we verify the L-condition for {X(t) : t G E} with respect 
to p and an appropriately chosen 9. That is, since the distribution functions 
Ft{-) have a density, they are continuous, and hence it suffices to show for 
an appropriate > there is a constant L < oo such that for every e > 

supP( sup \Ft{Xs) - Ft{Xt)\ > e^) < Le^ (72) 

t£E {s:seE,p{s,t)<e} 

We prove the L-condition holds assuming the density g of Z is uniformly 
bounded, and hence we have (j7Qp holding. The proof when g G -Zya(IR) is 
essentially the same, only the algebra changes, and hence the details are left 
to the reader. 

First we examine the situation when {y(t) : t > 0} is a martingale 
satisfying ([66]) with p = 1 and some /3 G (0, 1]. Applying (jTOj) to (j72]) we 
then have 

supP( sup \Ft{Xs) - Ft{Xt)\ > e^) < A, + B,, (73) 

teE {s:s(^E,p{s,t)<e} 

where 

e2 

= supP( sup \Xs-Xt\>—), 
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and 



Now 



£2 

= supP( sup \Xs - Xt\ > —). 

{s:se[{t-ee)yO,t]} 

Ae=SUpP( sup \Ys-Yt\> —), 

tGE 1 ^fc 

and hence Doob's martingale maximal inequality implies 

A, < snp2ke'^Ei\Y ^ - Yt\) < 2kCe-'+i (74) 

where the last inequality follows from ()66p with p = 1. We also have 

£2 

B,< sup P(\Y 1 -yt|>^) 

£2 

+ supP( sup lYs-Y^ 1 I > — ), 

tGE 1 t-eff VO' 

{s;s6[(t-e9)VO,t]} 

and using Markov's inequality, the martingale maximal inequality, and (66) 
with p = 1 as before, we have 

B,<8kCe-^+B. (75) 
Combining ([75]) . ([7ij) . and (f7S|) we have 

supP( sup \Ft{X,) - Ft{Xt)\ > e^) < lOA^Ce-^+l. (76) 

tGE {s:sGE,p(s,t)<t} 

Given our assumption that (j66|) holds with p = 1 and some /? G (0, 1], we 
take ^ = and hence (j75p implies we have the L-condition in (j72p with 
L = 10/cC < oo. 

Now we assume {Y{t) : t > 0} is a process with stationary independent 
increments satisfying (j66|) with p E (0, 1) and some /? E (0, 1]. Applying (j70|) 
to ([72]) we again have ([73]) . and as before 

£2 

A, = supP( sup |y,_y,|> ), (77) 

^ {s:sG[t,(t+e9)AT]} 

Since {Y{t) : t E -E} is a process with stationary independent increments 
and cadlag sample paths, an application of Montgomery-Smith's maximal 
inequality in jMS93| implies 

£2 

<3supP(|y 1 -Yt\>—). 

- tGE {t+eV)AT *' 20/c ^ 



26 



This maximal inequality is stated for sequences of i.i.d. random variables, 
but since {Y{t) : t G E} is a process with stationary independent increments 
and cadlag sample paths, for any integer n we can partition any subinterval I 
of E into 2" equal subintervals and apply |MS93] to the partial sums formed 
from increments over each of these subintervals. One can add auxiliary 
i.i.d. increments to form a sequence, but that is not necessary as at the 
n*'^each stage we need only work with the partial sums of the 2" increments 
of that stage. We then use [MS93| for an upper bound, and then pass via 
an increasing limit to what is needed, i.e. the desired upper bound is fixed, 
and hence is an upper bound for the limit. 

Thus by Markov's inequality, and our assumption of (I66p . we have 



A, < 3(20A:e-2)PsupE(|y i - Yt\P) < SCilOke'^fJ . (78) 

teE {t+e-<>)/\T 

We also have 

5, <supP(|y 1 -Yt\>—) 

+ supP( sup \Ys-Y 1 I > — r), 

{s:se[{t-e-5)\/0,t]} 

and using Montgomery-Smith's maximal inequality again we have 

£2 

B,<AsupP(\Y 1 -Yt\>—). 



Thus by Markov's inequality and (j66p 

B, < 4{40ke-^f supE{\Y i - Ytf) < 4C {40ke-^y . (79) 

teE {t+eS)AT 

Combining ([73]) , dTH]) , and ([79]) we have 

supP( sup \Ft{Xs) - Ft{Xt)\ > e2) < 7C7(40A:)f el'^f , 

teE {s:seE,p{s,t)<e} 

and hence the L-condition holds with L = 7C {40kY provided 6 = 2+2p • '-' 

Corollary 4. Let E = [0, T], and assume {Y{t) : t > 0} is a strictly stable 
process of index r S (0, 2] with stationary independent increments, cadlaq 
sample paths on [0, oo), and such that P{Y{0) = 0) = 1. Let X{t) = 
Z + Y{t),t > 0, where Z is a random variable independent of {Y{t) : t > 0} 



27 



and having density g{-) on M satisfying ()69p . Then, the empirical process 
built from i.i.d. copies of {X{t) : t G E} satisfies the CLT over C. Moreover, 
except for the degenerate cases when r = 1 and {Y{t) : t > 0} is pure drift, 
or Y{t) is degenerate at zero for all t €z E, the empirical CLT over C fails 
for these {Y{t) :te E}. 

Proof. The assertions about the CLT holding are immediate consequences 
of Proposition [3] once we check that {Y{t) : t > 0} satisfies (j66p . This follows 
from the comments immediately following (|66|) . and hence this part of the 
proof is established. 

To show that the CLT fails for the strictly stable stationary indepen- 
dent increment processes specified follows from an application of the Hewitt- 
Savage zero-one law, and the scaling property of such processes. The case 
r = 2 was previously established in jKKZlOj using a law of the iterated 
logarithm argument, which also applied to all fractional Brownian motions. 
Here we exploit the stationary independent increments of the processes to 
obtain a proof, and there are two other facts we need to emphasize at this 
point. The first is that {Y{t) : t > 0} strictly stable of index r £ (0,2] 
implies {S{t) = trY{j) : t > 0} is also strictly stable of index r. The second 
is that the process {Y{t) : t > 0} being non-degenerate and strictly stable, 
implies the random variables Y{t) and S{t) have probability densities for all 
t > 0. 

Now fix n > 1 and let Yi, • • • ,Yn be independent copies of Y. Let Q 
denote the rational numbers. Then, setting Sj{t) = t~Y{j),t > 0,j = 
1, • • • , n, we have 

P(card{Yi(t),-- - ,Yn(t)} = card{Si(t),--- ,Sn(t)}=n (80) 

for all t G Qn (0,oo)) = 1. 

Also, let Cq denote the countable subclass of C given by Cq = {Ct^y G C : 
t,y £ Q}. Then, as in the proof of Lemma 7 in [KKZlOj . to show the 
empirical CLT fails for {Y{t) : t G E} it suffices to show that 

p(A'^Q(yi,--- ,y„) = 2") = i, (81) 

where A^Q{Yi,--- ,y„) = card{C n {Yi, • • • ,Yn} : C G Cq}. Hence it 
suffices to show for every r,0 < r < n, and {ii, • • • , ir} ^ {1, • • • , n} that 

Pm„--- ,Yi^} e A^'^iYi,--- ,y„)) = L (82) 

Next define for every permutation tt = (ji, • ■ • ,jn) of {1, • • • , n} the event 

K = {w : Yj,{j,u;) < ■ ■ ■< Yj^{]-,u;) i.o. in k > 1}, (83) 
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and observe that 

Et, = {uj : Sj^{k,Lj) < ■■■ < Sj^{k,uj) i.o. in k > 1}. (84) 

Since there are only finitely many permutations and (|80p holds, P{Et^) > 
for some vr. Therefore, ()84p and the Hewitt-Savage zero-one law implies 
P{E-^) = 1. Moreover, since the processes li,--- ,l^n are i.i.d. it follows 
that P{Etj) = 1 for all permutations vr of {1, • • • , n}. 

Returning to (j82]) we take any permutation vr = (ji, • • • , jn) with ii = 
ji , • • • , ir = • Then for each uj and k such that (j83p holds we have a 
rational number q{uj, k) such that 

11 11 

Yjii]:^^) < ■■■< ^Jri-j:,^) < Qi^,k) < y,.^^^(-,L^) < ... < Yj„{-,uj), 

and hence 

{y,,(.,a;),... ,n,(.,^)} = Ci_^(,_^)n{yi(.,a;),... ,y„(-,a.)}. 
Since P{Et^) = 1, we therefore have ([52]) . which completes the proof. □ 

4 Applications to empirical quantile process CLTs 

The typical empirical quantile CLT of this section starts with a base process 
{Yt : t G i?}, and as before we define Xt = Yt + Z,t £ E, where Z is 
independent of {Yt : t G E} and Z has density g{-) on M. For the empirical 
process CLT's over C established in the previous section, we assumed g{-) 
was uniformly bounded on R, or in La(l^) for some a > 1. In order to prove 
our empirical quantile results, we assume a bit more about g{-), but these 
assumptions are not unusual, even for real-valued quantile CLT's. Moreover, 
keeping in mind possible application to a diverse collection of base processes, 
we have chosen to put the assumptions we require on g{-), but the reader 
should keep in mind that if the distributions ofYt,t S E, have densities with 
similar properties, then we could assume less about g{-). This is easily seen 
from the proofs, and basic facts about convolutions, and hence are left for 
the reader to implement should the occasion arise. 

Throughout we assume enough that the input process {Xf : t G E} 
satisfies the empirical CLT over C with centered Radon Gaussian limit on 
iooiExW) given by {G{t,x) : t £ E,x £ M}, where G{-, ■) is sample bounded 
on X M, and uniformly continuous with respect to its L2-distance there. 
Of course, as before a typical point {t,x) £ E xM has been identified with 
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Ct^x- Our empirical quantile CLT's in this setting will then be of two types, 
and in these results / will always be a closed subinterval of (0, 1). The first 
is that the quantile processes 



{V^(r^(t)-r,(t))/(t,r,(t)): n > 1} 



(85) 



satisfy the CLT in loa{E x /) with Gaussian limit process 



{G(t,Ta(t)): (t,a) E^x/} 



(86) 



and the second asserts that the quantile processes 



T^{t)): n> 1} 



(87) 



satisfy the CLT in ioo{E x /) with Gaussian limit process 



{ 



G{t,Ta{t)) _ 



{t,a) £ E X I 



} 



(88) 



Theorem 3. Assume that one of (i-iii) hold: 

(i) {Yt : t G E} is a centered sample continuous Gaussian process on a 
compact subset E of [0,T]'^ satisfying (f55]l . 

(ii) E = [0, r] and {Y{t) : i > 0} is a stochastic process with cadlag 
sample paths on [0, oo) such that P{Y{0) = 0) = 1. In addition, {Y{t) : t G 
E} is a martingale whose Li-increments are Lip-/3 for some /3 G (0, 1] , or a 
stationary independent increments process satisfying (I66p for some p G (0, 1) 
and /3 G (0,1]. 

(iii) E = [0, T] and {Yt : t G E} is a compound Poisson process built 
from the i.i.d random variables {1^ : A; > 1} having no mass at zero and 
Poisson process {N{t) : t >0} with parameter A G (0, oo) as in Proposition 



In addition, assume Xt = Yt + Z, where Z is independent of {Yt : t G E}, and 
Z has a strictly positive, uniformly bounded, uniformly continuous density 
function g on M. If {Yt : t G E} satisfies (i), (ii), or (iii), / is a closed 
subinterval of (0, 1), and we also assume that 



then the quantile processes of (|85]) and (f87|l built from the input process 
{Xt : t G satisfy the empirical quantile CLT with corresponding Gaussian 
limit as in (f86]l and dSSl) . 



3.2. 



fe^oo t^E 



lim snY>P{\Yt\ > b) = 0, 



(89) 
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Remark 3. It is easy to see at this point that the results of section three 
ahow us to apply Theorem 4.1 to obtain empirical quantile results of both 
types for fractional Brownian motions, the Brownian sheet, strictly stable 
stationary independent increment processes, martingales, and compound 
Poisson processes. The precise corollaries are easy to formulate, and hence 
are not included. 

Proof. If the base process {Yt : t £ E} satisfies (i), (ii), or (iii), and g is 
uniformly bounded on R, then Propositions 3.1-3.3 and Remark 3 following 
Proposition 3.2 imply that the resulting input process {Xf : t G E} satisfies 
the empirical CLT over C with centered Gaussian limit given by {G{t, x) : t € 
E,x £ M}, where G{-, •) is sample bounded on x R, uniformly continuous 
with respect to its L2-distance there, and has Radon support in iooiE x M). 
Furthermore, if Ht{x),t G E, is the distribution function of Yt, then Xt has 
probability density function 



Hence if g, the density of Z, is strictly positive, uniformly bounded, and uni- 
formly continuous on M, then it is easy to check that each of the densities 
f{t, ■),t S E, have the same properties. In particular, we have lim^^o ^'^PteE 
^^P\u-v\<s 1/(^7^) ~ fit,v)\ = 0, which is (f26]) . Hence Corollary 2.1 imme- 
diately implies the empirical quantile processes satisfy the quantile CLT's 
with Gaussian limit as indicated in ([86]) and ([88]) provided we show ([27]) 
holds. That is, it remains to verify that the densities f{t, ■),t E E, of the 
input process {Xt : t £ E} satisfy 



for every closed interval / in (0, 1). 

Now (|90p holds if we show that for any closed subinterval / of (0, 1) and 
all a > that 



First we show the left expression in ([9ip holds, so take a > 0. Then, for 




inf J{t,Ta{t))=ci>0 



(90) 



inf f(t,x) = Ca > and sup |TQ,(t)| < oo. 

teE,\x\<a teE,aeI 



(91) 
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every 6 > 

inf f{t,x)= inf / g{x — v)dHt{v) 

t&E,\x\<a t&E,\x\<aJ^ 

> inf / inf q(x — v)dHf(v) 

> inf qiu) inf / dHAv), 

- |«|<a+6^^ ' teE " 



and, since g satisfies (|89|) . there exists 60 > sufficiently large that 



inf / dHAv) > -. 
t^Ej-bo 2 



Therefore, we have 



inf f{t,x)>^ inf g{u) = Ca > 0. 

t£E,\x\<a 2 \u\<a+bo 



Now we turn to the second term in (j9ip . Since / is a closed interval of 
(0, 1) there is a 6^ G (0, i) such that I C (9,1-9) and 

sup P(|y(t) + Z\>a)< sup P (\Y{t)\ >^)+P (\Z\ > ^) 
ie[o,T] te[o,T] ^ ^' ^ ^' 

9 

where the second inequality follows from (j89p by taking a > sufficiently 
large. Hence for each t £ [0, T],a G / we have Ta{t) G [—«,«] and the right 
term of (|91|) holds. Thus (|9Up holds, and the theorem is proven. □ 



5 Additional quantile process CLTs with stable 
and Gaussian inputs 

Let E = [0, T], and assume {X{t) : t > 0} is a symmetric stable process 
of index r G (0, 2] with stationary independent increments, cadlaq sample 
paths on [0, 00), and such that P{X{0) = 0) = 1. Then, except for de- 
generate cases. Corollary |4] implies the empirical process built from i.i.d. 
copies of {X{t): t G E} fails the CLT over C, but it holds for i.i.d. copies 
of {X{t): a < t < T} over C^^^t] = {Ct,x : a < t < T,x £ M}, provided 
< a < T. Moreover, with E = [a,T] and I a closed subinterval of (0, 1), 
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Theorem [3] then imphes that the empirical quantile processes given in ()85p 
and (j87p . satisfy the CLT with limiting Gaussian processes as in (j86p and 
([88]) . respectively. 

As can be seen from the proofs, the difference in the results for X indexed 
by [0, T] versus [a, T] seems in large part due to the fact that the densities 
/(t, •) are not uniformly bounded on M as 1 1 0, and that X{Q) is degenerate 
at zero when t = 0. However, this is not the complete story, since in this 
section we will prove that the empirical quantile processes of (j87p with input 
process X on E = [Q,T] satisfy the CLT of ([88]) provided I is a closed 
subinterval of (0,1). As mentioned in the introduction, this extends the 
result of J. Swanson in ( [SwaOTj ) when I = {^} and in ( [Swall] ) for other 
fixed a G (0,1). 

To prove our result we need a number of lemmas. The first shows that 
if the input process is scalable, then certain information on the empirical 
quantile process on an interval, say, [1,2], yields information on an interval 
[0,5]. We phrase this in slightly more general terms, but ultimately it will 
be applied to quantile processes. 

Let {W{t,a) : t > 0,a G (0,1)} be a stochastic processes which is 
scalable in i.e. for some constant p G (0, oo) the process 

{W{ct, a) : t > 0, a G (0, 1)} and {c''W{t, a) : t > 0, a G (0, 1)} 

have the same law for all n > 1. Also assume P(VF(0,a) = 0) = 1 for 
all a G (0,1). Let Q denote the rational numbers, J = [1,2], A = [1 — 
a*, a*], ^ < Q* < 1, and for a subset i? of R we define Bq = i? n Q. 

Lemma 5. Let W be p-scalable. Fix < (5 G Q. For > 



E( sup 

Me(o,5]Q,aeAQ 



|T^(n, «)!"])< 



•E( sup \W{u,a)\'i]). 



(92) 



1 _ 2-P'i 



tie JqiOsAq 
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Proof. 



E sup \W{s,a)\'']=Esup sup \W{s,a)\'^ 

se(0,(5]Q,oeAQ j>l se(2-:'<5,2-(J-i)<5]Q,aeAQ 

oo 

< ^ E sup |VF(s,a)|« 

j^l se(2-J5,2-0-i)5]Q,aeAQ 
oo 

= ^ E sup \W{2-^5s,a)\'i 

se(i,2]Q,aeylQ 

oo 



Y^{2-^)P'iE sup |W^(s,a)|'^ 

j=l sG(l,2]Q,«eylo 

-E( sup \W{u,a)\i]). □ 



1 - 



Remark 4. Shortly we will apply this to the sequence of empirical quantile 
processes. That is, we apply this lemma to each of Wn{t, a) := y/n(^F~l {a) — 
F^"^(a)), where t > 0, a £ (0,1). Since the bounds obtained in Lemma [5] 
depend only on the scalability constant, c, all of our estimates will be uniform 
in n. 

We now prove that when the input process is a stationary, independent 
increment process with symmetric p-stable distribution (0 < p < 2), the 
empirical quantile processes uniformly satisfy the hypothesis of Lemma [5l 
For this purpose the next lemma is useful. 

Lemma 6. Let X be an arbitrary random variable. If Qai^) denotes any 
a-quantile for X, then —(/!_„(— X) is also an a-quantile for X. 

Proof. 

P{X > -q^_^{-X)) = P{-X < qi-a{-X)) > 1 - a 

and 

P{X < -qi-a{-X)) = P{-X > qi^a{-X)) > a. □ 

Theorem 4. Let {X{t): t > 0} be a symmetric r-stable process with sta- 
tionary, independent increments, and such that P{X{0) = 0) = 1. Then, the 
centered empirical quantile process built from i.i.d copies of {X{t) : t > 0} 
satisfies the hypothesis of Lemma O i.e., there exists a positive integer uq 
such that 

sup E[ sup y/ri\F~l{a) - Ff^{a)\] < oo, 

n>no t£jn,aeAQ 
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and hence for every e > there exists 6 > such that 

sup Pi sup \A^|Fj(a) - Fr'(«)l > e) < e- 

n>no te[0,5]Q,aeAQ 

Remark 5. From (llOOp . the no of the theorem can be taken to be 
no = mf{n > 1 : 2-('"Wi-"*)J-2)(^^^)2 < 

where AI! := , and Cr depends on the tail behavior of j as given 

I — a* 

prior to (f99|) . In addition, due to the scahng of our input process, Ff^{a) = 
t~F^^{a), and since a £ A imphes |Fj~"'^(a)| < F^^{a*) < oo, we have for 
< g < r that 

sup E[ sup V^\F-l{a)-Ff\a)\i]<oo, 

l<n<no tejQ.aeAQ 

i.e., maxi<„<„oSuptgj^^„g^^ |F-j^(a)| < maxi<j<„„ supo<j<r l^j WI < oo, 
and for < g < r the right term has a q^^ moment. Thus 

supE[ sup ^/n\F~l (a) - Ff^ {a)\''] < oo 
n>i teJQjoeAQ 

for < q < r, and the conclusion of Theorem [J] also holds for all n > 1 via 
an application of Lemma [5] with no = 1. 

Proof. First we note that the scaling property of the iid symmetric stable 
processes, {Xj{t)} immediately implies scalability with the same constant 
for all the processes 

{F-i(a) -F-i(a): t E [0,oo),a G (0,1)}. 

We'll obtain bounds on 

P{V^ sup \F-}{a)-FfHa)\>u). (93) 

strong enough to yield an no for which 

sup E[ sup ^/n\F~l{a) - Ff^{a)\] < oo, 

n>no tejQ,aeAQ 

At this point we can apply Lemma [5] to obtain the bound 

E( sup \F-Ua) - Fi\aW]) < (94) 
E( sup \F-}{a)-Fi\am. 
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An application of Chebyschev's inequality will then yield the Theorem. 

We break the proof into two parts. The first part covers the case when 
we have a lower bound on the densities of Xt for t G J = [1, 2]. In the second 
part we take care of the remaining case. Now for u > we have 

(i) 

sup (F„-;i(«) - Fr\<^)) > 

n 

= P{3t e JQ,a e AQ,y^ I u , >n(l-a)) 

X,(t)> — +F-\a) 

J— ^ Jn 



P{3t G JQ,a G AQ,Y,{lx,it)>u/V^+F,-He^) - ^(^(*) 

> u/V^ + Ff\a))) > n[(l - a) - P(X(t) > ^ + Ff^a))] 



n 

P{3t G jQ,a G AQ,Y,{lx,it)>u^+Ff\a) - Pi^it) 

> -^ + Fi\a))) > nP{Fi\a) < X{t) < -^ + F,-\a))). 



Also, since the density, /i, of ^(1), is symmetric about and unimodal 
( |Yam78] ). it is decreasing away from the origin. Hence, using < t < 
2, 

P{F^\a)<X(t)<^^F^\a))= / ^ /i(x)dx 



>( inf 

>h{F~\a*) + 



u , u 



Vi 

So, if < < C, since the density, /i, is decreasing away from the 
Jn 



origin, we have the inequality 



P{^ sup {F-}[a) - Fi\a)) > u 
1 " 

<n^\\Y.yx,it)>y-P{X{t)>y))\\j^ 
>/i(Ffi(a*) + C)n/2iA'). 



36 



Thus, for t £ Jq fixed, the continuity of P{X{t) > y) in y and the 
right continuity of Ixj{t)>y i^i U for ^ ^ j imphes 



1 " 

-^\\Y.i^x,ii)>y-p{m>y))\Wx^ 
1 " 

= ^11 ^(^w >y))ii^QX^^ 

with probabiHty one. 

Hence, if L> := l'^/"- MP-^ia*) + C), we have for < ^ < C, 

sup {F-}{a) - Ff\a)) > u) 
1 " 

< P(^|| Y.i^x,(t)>y - P{X{t) > y)) II j^xQ > ^n) (95) 

Since the summands, Ixj(t)>y " Pi^ii) > u) ^-re bounded by 1, and 
the CLT over C implies stochastic boundedness of the normalized norm 
in (195]) . we can use a result of Hoffman- J0rgensen, see pp. 164-5 of 
[HJ74] . to obtain for any q > 0, 

1 " 

■■= supE—\\Y,{lx,it)>y - P{X{t) > y))\\X,Q < oo. 

n V 

Therefore, for < n < C j \fn^ 

P{V^ sup {F-J{a) - F-i(a)) > n) < B^-^. (96) 

(ii) Now we deal with the case > C. In the computation below we 

don't use the particular form of the quantiles, F~l{a),Ff^{a), only 
the fact that they are quantiles. Note ([97|) below. Hence, by Lemma 
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El 

P{V^ sup \F~l{a)-FiHa)\>u) 

< P{V^ sup {F-J{a)-Ff\a))>u) 

+ P{V^ sup {Ff\a) - F-Jia)) > u) 

= P{V^ sup {F-J{a) - F,-\a)) > u) 

+ P{^/^ sup {-qt{l -a) + gn,t(l - «)) > u). (97) 

Thus, the second term can be treated the same as the first term. For 
the first term we have 

P{V^ sup {F-}{a) - Ff\a)) > u) 

= P{3t G Jq, a G Aq, F-Jia) > ^ + F^-^a)) 



P(3t G Jq, a G Aq, 31, #/ = [n(l - q)J , Xj{t) 



>^ + Fr^(a),VjG/) 

<P(3teJQ,a€AQ,3I,#I= [n{l-a*)\,Xj{t) 
> ^+Fi\a),yj el) 



in 

and again by Lemma El since F^^{a) > Ff^{l - a*) for all a G Aq, 
l)p(3t G JQ,X,it) > - Ff\a*), 



[n{l - Q*)J 
j = l,...,Ln(l-a*)J) 

n 

[n(l — a 

£ Ir^pumij, > ^ - -Ff'K))] (98) 

Now, for our stable process it is known, see Proposition 5.6 of |LT91] . 
that there exists a constant, Cr, such that 

P{\\X\\j^ >V)< CrV-\ 
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Hence, by mh, if ^ > C > 2F.-\a*), we have for A: := - 
Vn 1 

P{V^ sup {F-l{a)-F-\a)) >u) 

< r g ^ f V^Nrl Ln(l-a')J ^ rW^.r-Ln{l-Q*)J .^g. 

- ^1 - a* ''"^ n J ^ u ^ ' ^ ^ 

Therefore, taking n/^H >C = 2\rV 2Ff ^(a*), ([Ml) imphes 

P{V^ sup (F-i(a)-Fr^(a))>^)<[^]'^Ln(i-*)J 



and hence n sufficiently large implies 

sup (F-i(a) - Fr'(a)) > n) < u-\ (101) 

teJoiQeAo 



Since the same estimates apply to the second term in (j97j) , we have by 
putting the two parts together that 

E[ sup V^\F-l{a) - Ff\a)\] 

POO 

<2 sup {F-^{a) - F^^{a)) > u)du 

Jo teJQ,a£AQ 



< 2[1 + / sup {F-J{a) - F-i(a)) > u) du 

JO tejQ,a£AQ 

POO 

+ PiV^ sup {F-Jia)-Fi\a)) >u)du] 

Bo 1 [°° I 

< 2[1 + / -^du+ ^ dn] < oo, (102) 



provided n is sufficiently large, B2 and D are as in ()96p . and C = 
2A.f. V 2F^^{q*). Thus the hypotheses in the LemmaOare uniformly 
satisfied. Hence the theorem is proved. □ 
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Our next lemma is important in that it allows us to switch back and 
forth between supremums over countable and uncountable parameter sets. 
For example, one consequence is that the conclusion of Theorem U] can be 
strengthened to hold for all t G J and a £ A provided we ask that the 
processes {Xj{t) : t > 0},j > 1, have cadlag sample paths with probability 
one. The lemma is as follows. The sets A and Aq are as above, but [0,T](Q 
also includes the point T, even if it is irrational. 

Lemma 7. Let {X{t) : t > 0} be a symmetric r-stable process with station- 
ary, independent increments, cadlag sample paths, and such that P{X{0) = 
0) = 1. Let < T < oo, Q the rational numbers, and define [0,T](Q = 
([0, T] n Q) U {T}. Then, the empirical quantile process r"(i) built from i.i.d 
copies of {X{t) : t > 0} with cadlag paths on a complete probability space 
has right continuous paths on [0, T] with probability one, and is such that 

Pi sup |r^(t) -r„(t)| = sup |T^(t)) - r«(t)|) = L (103) 
te[0,T],aeA tG[0,T]Q,aeAQ 

Moreover, for each t S [0, T] and n > 1 with probability one the empirical 
quantile process T^{t) is left continuous in a E (0, 1). 

Proof. First we observe that for t > the distribution function F{t, x) has 
strictly positive density given by 

f{t,x) = {2-k)~^ / exp {-ct\u\''} cos{xu)du, c > 0, (104) 
Jr 

and hence F{t, x) is strictly increasing and continuous in x G M. Thus 
Ta{t) = Ff (q^) is continuous in a G (0,1), and by its definition, we also 
have T^{t) = F~1{q) left continuous in a G (0,1). In particular, the claim 
following ()103p holds. Moreover, for every a G (0,1), P(tq,(0) = r"(0) = 
0) = 1 since we are assuming P(Xj(0) = 0) = l,j > 1. Therefore, we have 

P{ sup \T^{t) - T^{t)\ = sup |r^(t))-T„(t)|) = l. (105) 

te[0,T],aeA te[0,T],aGAQ 

We also have Ta{-) continuous in t on [0, cxd), since scaling easily implies 
T^(t) = t-^Ta{l) for all t > 0. Thus (fTMl) follows from (fT05]l provided we 
show T^(t) is right continuous on [0,T) with probability one, i.e. we then 
would have 

P{ sup \T^{t) - T,,{t)\ = sup Kit)) - T^{t)\) = I. (106) 

tG[0,T],aeAQ te[0,T]Q,aeAQ 
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To verify the right continuity of r2'(-)) and (|106p . we use the right conti- 
nuity of the paths of the processes Xi, • • • , X^- We do this through the fol- 
lowing observation. That is, given real numbers xi, ■ ■ ■ , x„, let xi^i, • • • , xi^„ 
be an ordering of these numbers such that xi^i < • • • < xi^n- In case there 
are no ties among {xi, ■ ■ ■ , x„} this ordering is unique, and when there are 
ties, we choose the ordering based on the priority of the original index among 
the tied numbers. We then refer to xi^i < • • • < as the order statistics 
of {xi, • • • , Xn}- Now take pairs {xi,yi), • • • , {xn,yn) of real numbers such 
that supi<j<„ \xj — Ujl < 6. These are called the initial pairs of the two sets 
of n numbers. We will now verify by induction that the corresponding order 
statistics formed from these sets also satisfy 

sup \xij - yij\ < 5. (107) 

The case n = 1 is obvious, so assume the result holds for all sets with 
cardinality less than or equal to n — 1. Then, assume that in the initial 
pairings, xi^i and yi^i are paired with, say yi^k and xi^i, respectively. Hence 
we have < < yi^k, \xi,i - yi,k\ < ^ and - xi^i\ < 5. If 

xi,i < yi,i, then from the above we have xi^i < yi^i < yi^k-, and hence 
|xi^i — yi^k\ < ^ implies \xi^i — yi^i\ < 5. Similarly, if xi^i > yi^i, then from 
the above we have yi^i < xi^i < xi^i, and hence \yi^i — xi^i\ < 5 implies 
l^^i,! — yi,i\ < ^- We also have \xi^i — yi^^l < ^- That is, if xi^i < yi^k-, then 
we have xi^i < xi ^ < yi k and hence |xi^i — yi^k\ ^ ^ implies |xi^; — yi^k\ < 
5. Similarly, if xi^i > yi^k-, then we have yi^i < yi^k < xi^i and hence 
\yi,i - xi^il < 5 implies |xi,/ - yi^k\ < S. 

To finish the proof of ()107p we apply the induction hypothesis to the 
set of n — 1 pairs determining xi_2 < • • • < xi^„ and yi_2 < • • • < yi,n, with 
ixi,i,yi,k) being a possibly new pair, and the remaining n — 2 pairs are those 
originally given. Note that the induction hypothesis applies to these pairs, 
since we have shown |xi^/ — yi^k\ < Thus (I107P holds. 

To verify the right continuity of t^(-), and hence that (jl06p holds, we 
note that since the i.i.d. processes Xi,--- ,Xn are cadlag on [0,oo) with 
probability one, there is a set r^i C Q such that P{i^i) = 1 and for every 
t G [0, T), e > 0, there is a (5 = t, e, n) > such that uj £ implies 

sup \Xj{s) - Xj{t)\ < e. 

i<j<n,t<s<{t+S)/\T 

Therefore, (|107p implies the order statistics Xi^i{s) < ••• < Xi^„(s) and 
-^^1,1(0 < • • • < ^i,n(t) obtained from {Xi{s), ■ ■ ■ ,X„(s)} and {Xi{t), ■■■ , 
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Xn{t)} are such that 

sup \Xij{s) - Xij{t)\ < e. (108) 

^<j<n,t<s<{t+S)AT 

Since e > is arbitrary, we thus have that the order statistic processes 
{Xij{t) : t G [0, r)},j = I,-- - ,n, are right continuous on ili, and hence 
with probabihty one. 

Now for < a < l,n > l,t G [0,oo) we have T^{t) = inf{x : Fn{t,x) > 
a}, and hence 

where j{a) = min{A: : 1 < k < n,k/n > a} is independent of t £ E. Thus 
for ah uj £ rii we have T^(t) right continuous in t G [0, T). Hence the lemma 
is proven. □ 

Corollary 5. Let {X(t): t > 0} be a symmetric r-stable process with 
stationary, independent increments, cadlag sample paths, and such that 
P{X{0) = 0) = 1. Also, assume the empirical quantile processes T^{t) are 
built from i.i.d copies of {X{t) : t > 0} with cadlag paths. Then, there 
exists an integer uq > 1 such that for every e > there is a (5 = 6{e) > 
satisfying 

sup P( sup V^Kit) - Ta{t)\ > e) < e. (109) 

n>no te[0,(5],oeA 

Proof. The proof of the corollary follows immediately from Theorem |4] and 
Lemma [71 □ 

The next theorem shows an empirical quantile CLT holds on [0, T] for the 
symmetric stable processes discussed here, which contrasts with the remarks 
at the end of section 3 showing that the empirical CLT for such processes 
fails. It also extends the results of |Swa07j and [Swallj for Browning motion, 
showing the CLT is uniform in a G /, where I is a closed subinterval of (0, 1). 

Theorem 5. Let {X(t): t > 0} be a symmetric r-stable process with 
stationary, independent increments, cadlag sample paths, and such that 
P{X{0) = 0) = 1. Also, assume the empirical quantile processes T^{t) are 
built from i.i.d copies of {X(t) : t > 0} with cadlag paths, and / is a closed 
subinterval of (0, 1). Then, the quantile processes 

{V^(r:(t)-r„(i)): n > 1} (110) 

satisfy the CLT in -^c>o([0,X'] X with. cGiitered Gcnissiajii limit process 

{W{t,a): {t,a) G [0,T] x /} , (111) 
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where W{0, a) = 0, a G I, 

W{t, a) = ^fll'^ill] , {t, a) e (0, T] x /, (112) 
and for (s,/3), {t,a) G (0, T] x / the covariance function is given by 

f{s,rfi{s))f{t,Ta{t)) 

Remark 6. Since the process {X{t) : t > 0} in Theorem [5] is scalable with 
parameter i, it is easy to check that Ta{t) = t~rQ,(l) for {t,a) € [0, cxd) x 
(0, 1). In addition, since the density of X{t) is strictly positive for each 
t > 0, it is easy to check for t > 0, x £ M that 

f{t,x)=t-^f{l,Xt-'r). 

Thus ioTt>0,ae (0,1), 

f{t,T^{t))=t-"-f{l,T^{l)), 

and for {t,a), (s,/3) £ {0,T] x /, (fTT3D becomes 

j(l,-^/3(l))/(l,Ta(l)) 

Furthermore, since 1^(0, a) = 0, a G /, we also have (jll4p when (f, a), (s, /3) 
e [0, T] X I. 

To get the covariance for the limiting median process, set / = {i}. Then, 
since Ti{t) = for all t G [0,T], we have the limiting Gaussian process such 

2 

that P{W{0, |) = 0) = 1 and for s,t E [0,T] its covariance is 

. . l^ . l^^ 1 1 P(X (s) < 0, X (t) < 0) - ^ , , 

The density f{t,x) is as in ()104p . and hence 

/(t,0) = i^^£t™.,>0, (116) 

2lT{ct)r 

which implies for s,t £ [0, T] that 

E(^(., l)wit, i)) ^ [^(^(^) < o, xit) < o - 1]. 

(117) 
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In these examples the sample paths of the input process X are assumed to 
be cadlag, and when r = 2 they could be assumed to be continuous as X 
is then a Brownian motion. Hence, after the proof of Theorem [5] we will 
examine the consequences of these path properties for the quantile CLT. 
At this time we also will discuss results comparable to those for sample 
continuous Brownian motion when the input data comes from any sample 
continuous fractional Brownian motion. 

Proof. Let 

Wnit,a) = V^{{T^{t)-To,{t)),t G [0,T],a G /,n > 1. 

Then, P(Wn{0,a) = 0) =1 for a £ I,n > 1, and the finite dimensional 
distributions of Wn converge to the centered Gaussian distributions given by 
the covariance function in (jll3p for t € {0,T],a G I. Hence Theorem 1.5.4 
and Theorem 1.5.6 of |vdVW96] combine to imply the quantile processes 
{Wn:n>l} satisfy the CLT in ^qq([0,X'] X /), where the limiting centered 
Gaussian process has the covariance in (|113p . provided for every e > 0,r] > 
there is a partition 

[0,T]x I = utiEi (118) 

such that 

limsupP*( sup sup \Wn{t,a) - Wn{s,P)\ > e) < ry. (119) 

n^oo l<i<k (t,a),(s,l3)<^E, 

Since / is a closed subinterval of (0, 1), there is an a* G (^, 1) such that 
/ C ^ = [1 - a*, a*]. For (5 > and Ei = [0, (5] x / observe that 

P*{ sup \Wn{t,a)-Wn{s,p)\>'^) 
< 2P*( sup \Wn{s,a)\>'l^). 

s&[0,S],a&I ^ 

Hence, and (fT09D imply there is a 5 = 5{^) such that 

limsupP*( sup \Wn{s,a) - Wnit, > ^) < 2(^) < ^. 

n-s-oo {s,a),{t,l3)eEi ^ 4 2 

(120) 

Now Theorem [3] above implies the CLT for {Wnii-, ol) : (t, a) G [(^, T] x /} 
in lcyQi\b^T\ X /), and hence Theorem 1.5.4 of |vdVW96] implies that there 
is a partition [(5, T] x / = U^^2-^j such that 

limsupP*( sup sup |VF„(s,a) - VF,(t,/3)| > ^) < (121) 

n-!>oo 2<i<fc (s,Q),(t,/3)G£;i ^ ^ 
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Combining (|120p and ()12ip we have (jllSp and ()119p . and hence the theorem 
is proved. □ 



Now we turn to the question of how special sample path properties of the 
input process {Xt : t £ E} influence our quantile CLT's. To be more spe- 
cific, recall that the CLT results we have established for empirical quantile 
processes, hold uniformly in the space looiE x /), and the limiting Gaus- 
sian process {W{t,a) : {t,a) G E x /}, almost surely, has a version with 
paths which are bounded and uniformly continuous with respect to its own 
L2 distance dw on E x I. In particular, this guarantees that the measure 
induced by the Gaussian process on loa{E x I) is supported on the subspace 
Cl^{E X I) of £00 {E X I), where the subscript L2 is written to indicate the 
topology on £' X / is that given by the Gaussian process L2 distance. Hence, 
if the input process {Xt : t £ E} is assumed sample continuous on {E,ei), 
where ei is a metric on E, when does our quantile CLT with a £ (0, 1) fixed 
hold on the space of ei continuous paths? li E = [0,T] with metric the 
usual Euclidean distance, and the input process has cadlag sample paths on 
[0,T], a similar question can be asked if the quantile CLT holds in some 
related space of functions. Since processes with continuous paths or cadlag 
paths are typical of many examples throughout probability and statistics, 
these are natural questions, but they also relate to some recent results of 
Jason Swanson. That is, he established a CLT in the space of continuous 
functions on [0, T] for the median process obtained from sample continu- 
ous Brownian motions in |Swa07] . and for other individual quantile levels 
a G (0, 1) in [Swallj . These results will follow from our next theorem, and 
are established in a remark following its proof. 

Since the empirical quantile processes have jumps as a ranges over (0, 1), 
to state our theorem providing some facts related to these questions, we need 
the following function spaces. If ei is a metric on E we set 

Cei{E) = {z : z is continuous on (E, ei)}, (122) 

if -E = [0, T] we assume ei is the usual Euclidean distance and let 

Bi([0,r]) = {z: zis cadlag on [0,T]}, (123) 

where right and left limits are taken with respect to ei on [0, T], and for 
/ = [a, b] a closed subinterval of (0, 1) we set 

]D)2(/) = {z : z is left continuous on (a, b], and has right limits on [a, b)}, 

(124) 
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where right and left hmits are taken with respect to the usual Euclidean 
distance 62 on /. We also define the closed subspaces of ^00 {E x /) given by 

Ce,{E)®B2{I) = {f{;a) G CeAE) E I and /(i, ■) G Vt G [0,r]}, 

(125) 

and 

Bi{[Q,T])^B2{I) = {/(-,«) G Bi([0,r]) Va G 7 (126) 
and/(t,-)eB2(/) VtG [0,r]}. 

Both Ce^E) (g) B2{I) and ]D)i([0,r]) (g) ^2(1) are closed subspaces of 
£^{E X /)). 

Theorem 6. Let {Xt : t G be the input process for the empirical quantile 
processes defined for t e E,a e (0, 1), n > 1, by 

Wnit,a)=V7iiT^it)-Ta{t)), 

and assume they satisfy the empirical quantile CLT in i^o {E x I) with Gaus- 
sian limit {W{t,a) : {t,a) G E x I}. Let dw denotes the L2 distance of W 
on E X I, and assume the identity map j on £■ x / is continuous from the 
ei X 62 topology on £^ x / to the dw topology, and that Ta{-) G Ce^iE) for 
every a G /. Then, we have: 

(i) If {Xt iteE} has version with paths in Cg^iE), then the empirical 
quantile CLT holds in the Banach subspace Cei{E) (g) ID2(7) of ioo{E x I). 
In particular, if a G (0, 1) is fixed, then the CLT will hold in the space of 
continuous functions Ce^ (E) with the topology that given by the sup-norm. 

(ii) ltE= [0, T] and {Xt : t e E} has version with paths in ]D)i([0, T]), we 
have the empirical quantile CLT holding in the Banach subspace Bi([0, T])(8) 
]D)2(/) of £oo([0,T] X I). Hence, if a G (0, 1) is fixed, then the CLT will hold 
in the space of functions ]D)i([0,r]) with the topology that given by the 
sup-norm. 

Proof. If {Xt -.teE} has a version with paths in Ce^iE), then taking i.i.d. 
copies of this continuous version to build the quantile process, the proof of 
Lemma 5.3 implies one has with probability one that Wnit, a) is continuous 
on (E, ei) for each a G (0, 1). In addition, for each n > 1, by Lemma 5.3 we 
have a — >■ T^{t) is in B2(/) for all t e E, and therefore 

P{Wn{;-)^CeAE)(^B2{I)) = l. (127) 

Moreover, since we are assuming the identity map j is continuous from the 
ei X 62 topology to the dw topology on E x I, with W having a version 
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with paths in Cl2{E x /), it foUows from the fact that the composition of 
continuous maps is continuous that W has a version such that 

eCe,(^)»P2(/)) = 1. (128) 

Combining (jl27p and (jl28p , an easy appUcation of the portmanteau theorem 
imphes (see Theorem 1.3.10 of [vdVW96j for the details) the CLT will hold 
on Cei(-E') ^ 1152 (/) with the topology that is given by the sup-norm. Thus 
(i) holds. 

The proof of (ii) is entirely similar, since the assumptions of (ii) and 
Lemma[f|imply that (fT271) holds with Ce^{E)0B2{I) replaced by Di([0, r])(g) 
D2(-^)- Moreover, since we always have 

Cei([0, T]) ©2(1) C Bi{[0, T]) ® D2(/), 

and the argument for (jl28p is valid under (ii) , we have that (I127P holds with 
this replacement. Hence (ii) is verified as before. □ 

Remark 7. In this remark we provide specific applications of Theorem[6l Our 
first application assumes the input process X is a cadlag symmetric r-stable 
process with stationary independent increments on [0,T] with P{X{0) = 1, 
and shows that under these conditions the empirical quantile CLT holds 
in the Banach space Bi([0, T]) with the sup-norm. The special case r=2 
implies X is Brownian motion, and if the quantile CLT is built from i.i.d 
sample continuous Brownian motions, then we will also see that for fixed 
a G (0, 1) the empirical quantile CLT holds in the Banach space Cei([0,T]) 
with the sup-norm. As we mentioned earlier, this implies the quantile CLT 
for medians in |Swa07) . and for other individual quantile levels a G (0,1) 
in [Swallj . A major step in these results will be the use of Theorem [5] to 
establish the empirical CLT. 

In the second application the input process X is a sample continuous 
fractional Brownian motion, and here for fixed a G (0, 1) we again have the 
empirical quantile CLT in Cei([0,T]) with sup-norm. However, for this class 
of examples we will only outline the necessary arguments as they are much 
the same as those for the stable processes. Hence we now turn to that case. 

First we observe that if X is a cadlag symmetric r-stable process with 
stationary independent increments on [0,cxd] with P(X{0) = 1, then for 
T > fixed and / a closed subinterval of (0, 1) we have the empirical quantile 
CLT of Theorem [5j Furthermore, since {X^ '■ t > 0} is equal in distribution 
to cr{Xt: t > 0} for c > 0, it easily follows that Ta{t) = trTa{l) is jointly 
continuous in {t,a) £ [0, oo) x (0,1). Hence Theorem [U] implies the claims 
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made above for the stable process inputs provided we show the identity map 
j on [0, T] X / is continuous from the EucUdean topology to the dw topology 
on [0, T] X /, where ()114p and (s, (t, a) e [0, T]xl imphes the L2 distance 
dw is given by 



,9 // / sr(B — B) tr(a — a 



_2 i.i [P{Xs<Ms),Xt<To,{t))-af3] 

Hence, fix {t,a) G [0, T] x I, and assume e((s, /?), (t, a)) = ei(s,t) + 
62(0:, /3) — )• 0. Then, if t = it is obvious that the identity map j is 
continuous at (0,a),a G /, as asserted, i.e. e((s, /3), (0, a)) — )• implies 
dvK(('S,/3), (0,a)) — )• since the density is strictly positive and con- 

tinuous in X G M and t^(1) — )• Ta(l) as /? — )• a > 0. Moreover, for 
t G (0, T] and a £ I fixed, the continuity is again obvious provided we 
show e((s, (3), (t, a)) — implies 

PiX, < T^is),Xt < r„(t)) ^ a. (129) 

To verify (|129p we observe that if An and i?ri, are sequences of ran- 
dom variables on the same probability space with plimj^_j,oQAn = A and 
plimjj_j.oQBn = B, then the random vector {An,Bn) converges to {A,B) in 
probability, and P{An < x, Bn < y) converges uniformly to P[A < x,B < y) 
for (x, y) G provided [A, B) has a continuous distribution function. In 
particular, ii A = B, then (j4„,i?„) converges in probability to {A, A), and 
P{An < x, Bn < y) converges uniformly to P{A < x,B < y) = P{A < xAy) 
for {x,y) G provided A has a continuous distribution function. Now Xt is 
continuous in probability on [0,T], and t^(s) = s~r^(l) is jointly continuous 
in s G [0, T] and /3 G (0, 1), so setting An = A = Xt with t > 0, a G / fixed, 
and Bn = Xs„ , Sn ^ t, f3n ^ a, B = Xt, we then have 

lim P{Xs„ < T^„(Sn),Xt < r„(t)) = P{Xt < Ta{t)) = a. 

Since the sequences s„ and /3„ with the stated properties are arbitrary, we 
have (I129p . Thus the claims regarding the stable processes of this remark 
are established. 

Now we turn to the application of Theorem [6] to fractional Brownian 
motions. Hence let {Xt : t > 0} be a centered sample continuous 7-fractional 
Brownian motion for < 7 < 1 such that Xq = with probability one 
and E(X|) = t^"^ for t > 0. Take / a closed subinterval of (0,1), and 
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assume the empirical quantile processes built from i.i.d. copies of X with 
continuous paths are given for t > 0,a £ (0,1), n > 1, as in (jllOp . Since 
{Xct '■ t > 0} is equal in distribution to c^{Xt : t > 0} for c > 0, it follows 
that Ta{t) = f^Tail) is jointly continuous in (i, a) £ [0, oo) x (0, 1). Moreover, 
the analogue of the argument given above for stable processes implies the 
identity map j on [0, T] x / is continuous from the Euclidean topology to 
the L2 distance dw on [0, T] x /, where in this situation 

dl,{is,/3),{t,a)) = 27r{.27exp{r|(l)}(/3-/32)+t27exp{r2(l)}(a-a2)_ 

(130) 

2.^t^exp{l(r2(l) + r|(l))}[P(X, < T(s{s),Xt < r«(t))-a/3]} 

for (s, P),{t,a) G [0, T] x /. Hence we will have empirical CLT results 
as in part (i) of Theorem [6] with E = [0,T], provided we can show we 
that Theorem [5] applies when the input process X is a fractional Brownian 
motion. 

Except for j = ^, X does not have independent increments, but the 
results of Theorem H] still hold in this setting. Moreover, from the proof of 
Lemma [7] the empirical processes are such that 

PiWni;-) (^CeAE)(E)B2{I)) = l. (131) 

Furthermore, Corollary [5] and Theorem [5] hold when the input processes 
are sample continuous fractional Brownian motions, and hence the proof 
of Theorem [5] implies the empirical quantile CLT in £oo([0,T] x I) in this 
setting. The limiting Gaussian process W{t, a) has covariance as in ()113|) 
for s,t £ (0, T] and zero for s or t equal to zero, and hence the covariance for 
W is as given in (I130p . Hence the empirical quantile CLT results in part (i) 
of Theorem [6] hold as indicated with E = [0, T] and a fractional Brownian 
motion as the input process. 

In particular, if a G (0, 1) is fixed, then the CLT would hold in the space 
of continuous functions on [0, T] with the topology that given by the sup- 
norm. The limiting Gaussian process then has covariance as in (jl30p with 
(3 = a. Hence if a = /3 = ^ and s,t £ (0,T], the covariance of the limiting 
Gaussian process is 

wm/r ^u/r. PiXis)<0,X{t)<0)-l . _ EiXsXt). 

E(W{s, -)W[t, -)) = — 7 = s't^ sm (— -) 

^ ^ 2' ^ ' 2" f{s,0)f{t,0) ^ sTH 

where 2K{XsXt) = s^"' + t^'^ — k — tp'^, and the second equality follows 
from a standard Gaussian identity. When 7 = 1/2 and the input process is 
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standard Brownian motion, this gives the covariance in [Swa07] . Thus the 
claims of this remark are estabhshed. 

Remark 8. The Brownian sheet also has a scaling property, and continuous 
paths. Hence by using the methods of the previous remark, results of the 
type discussed there should also hold for the sheet. We have checked these 
results when the input process is the 2-parameter sheet on [0,r] x [0,r], 
but the details differ very little from what is done in the previous remark, 
so they are not included. 
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